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Every Bose-Einstein condensate is in a highly entangled 
state, as a consequence of the fact that the particles in a 
condensate are distributed over space in a coherent way. It 
is proved that any two regions within a condensate of finite 
particle number are entangled. This entanglement does not 
depend on the distance between the two regions. Criteria 
for the presence of entanglement are derived in the context 
of interference experiments. For separable states there is a 
trade-off between fluctuations in particle number and inter- 
ference visibility. 

Entanglement is one of the dramatic non-classical fea- 
tures of quantum physics. There are states of composite 
quantum systems that cannot be decomposed into prob- 
abilistic combinations of product states; there is no way 
of writing the density matrix of such a state p in the form 
p = ^pia\® a\® <jI^, where a\ is a state of the k-ih 

i 

subsystem and the pi > Q are probabilities. Such states p 
are called inseparable or entangled. The properties of an 
entangled system are not completely determined by the 
properties of its parts. Only entangled states can exhibit 
quantum non- locality 0|. In this case there is no way 
of reproducing the predictions of quantum physics with 
classical systems, unless there is instantaneous commu- 
nication over arbitrary distances. 

In the last few years, entanglement has been stud- 
ied extensively in the context of quantum computation 
and quantum communication |^]. Methods for creat- 
ing entangled states have been found in various physi- 
cal systems, including non-linear optics, ion traps, cav- 
ity quantum electrodynamics and nuclear magnetic reso- 
nance 1^. A fascinating open problem is whether entan- 
glement plays a significant role in natural physical sys- 
tems. First steps towards addressing this question were 
made in the study of interacting spin systems Q . The ef- 
fect of decoherence on large-scale spin entanglement has 
also recently been investigated in this context . 

Bose-Einstein condensation is another genuine quan- 
tum phenomenon which has recently attracted a large 
amount of experimental and theoretical attention In 
a Bose-Einstein condensate (BEC), all particles (which 
have to be bosons) are in the same quantum state. 

Here I would like to draw attention to the fact that by 
its nature every BEC is a highly entangled state. Even 
when there is no entanglement of the internal states [Q, 
there is always entanglement in the external degrees of 
freedom. The discussion will be focussed on the entan- 
glement between different spatial modes, but analogous 
statements can be made for momentum modes. It will be 



shown that any two regions inside a BEC of finite particle 
number are entangled with each other. The detection of 
this entanglement via interference will also be discussed. 

It is convenient to use the language of non-relativistic 
quantum field theory. The particles composing the 
BEC are described by a bosonic quantum field operator 
'(/'(x) satisfying the commutation relation [ip{x), — 
5^'^\x — y)!, where x and y are points in three- 
dimensional space. The time dependence of il>{x) is not 
essential for the present purpose. Throughout the paper 
operators will be distinguished by hats. 

A BEC of N particles is described by the state I^E") = 
:^7^(/^)^|0)j where |0) is the vacuum state and the cre- 
ation operator y can be written in terms of '(^{x)^ as 
P = J (fix f* {x)i()^ {x). The mode function f{x) is nor- 
malized: J (fix \f{x)\'^ = 1, where the integrals are over 
all of space. It is easy to check that / fulfills the commu- 
tation relation [/, f^ = 1. 

Let us consider a certain volume inside the BEC and 
study its entanglement with the rest. Denoting the vol- 
ume by A, one can define mode operators 

a = y d^xf{x)4'{x) 

xeA 



b=^ J d\f{x)i:{x), 

xeA 



(1) 



where A is the complement oi A, p = J (fix\f{x)\'' 
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and q ^ 1 — p = J (fix\f{x)\'^ . Then one has / 
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y/pa+^b, and a and b fulfill the following commutation 
relations: [a,b] = [a,b^ = 0, [a, a^] = [b,b^ = t. 
The state of the BEC can now be written as 
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(2) 



which clearly shows the entanglement between spatial 
modes a and b. These modes can have different physical 
meanings. For example, a could be the interior of the 
BEC and b the rest of the universe including the BEC's 
boundary region; a and b could also each contain one half 
of the BEC. 

The nature of the entanglement in (^ is easy to see. 
The N particles can be distributed in various ways over 
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the two modes a and h. The state of the BEC is a coher- 
ent superposition of all these possibilities. This clearly 
generalizes to more than two modes. Note that the par- 
ticle numbers in d and b have a binomial distribution. 

The present kind of entanglement cannot be detected 
by local measurements, i.e. measurements acting on 
modes d or h separately. If one is restricted to local mea- 
surements in the basis of particle number, then is 
indistinguishable from a separable mixture of the various 
terms without any fixed phase relation between them. 
One would therefore require local measurements in bases 
of states corresponding to superpositions of different par- 
ticle numbers, which is impossible. However, this does 
not mean that the entanglement in (H) is unobservable. 
It can be detected by joint measurements on the two 
modes, for example by interference experiments. This 
point will be discussed in more detail below. 

Let us now study the entanglement properties of BECs 
in more detail. First suppose that instead of a Fock state 
with exactly N particles the BEC is in a coherent state 
of the form Q 

This can formally be rewritten as 



la) =e-HV2 



|0), 



(4) 



where the product is over all points in space. The state 
is clearly a product state with respect to the point- 
like spatial quasi- modes ^{x), since the vacuum does not 
contain any entanglement |^ . 

As a consequence all states that can be written as con- 
vex combinations of coherent states are also unentan- 
glcd. In particular this applies to a Poisson distribution 
of number states \n) = '■■^-1^. |0): 
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|n)(n|. 
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which can be rewritten as a mixture of coherent states 
with fixed amplitude \/\ but random phase 
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pp = j d(j}\V\e"^){VXe"l'\ 




(6) 



as can easily be shown by expanding the coherent states 
in the number state basis: \a) = e"'"' ^T^l"-) ^'^'^ 
performing the integration over (j). 

Note that a state with a thermal distribution of par- 
ticle numbers is also not entangled, because it can be 
written as a convex combination of coherent states with 
a probability distribution that is a Gaussian in the co- 
herent state amplitude: 



PT = {l-t)^e\n){n\= /d^ae-'^H (7) 



where the integration is over the complex plane and 
/3 = (1/i) — 1. Again the identity can be shown by ex- 
pressing the coherent states in the number state basis 
and performing the integration. 

The above results suggest that the fixedness of the par- 
ticle number plays an important role in ensuring that the 



state 



■|0) is entangled. This leads naturally to the 
question of subsystems. If one considers a volume con- 
taining just part of a BEC, then the particle number in 
this volume is no longer fixed. Does the reduced state 
corresponding to such a volume still exhibit spatial en- 
tanglement of the above type? It will now be shown that 
the answer is yes, and that the entanglement disappears 
only in the limit of a finite subsystem of an infinitely 
large BEC. 

Let us partition space into three regions. A, B and 
C. For concreteness, imagine A and B to be well 
within the region occupied by the BEC. A and B may 
be adjacent or separated. We will trace over the re- 
gion C and study the entanglement between the regions 
A and B. Defining modes d, b and c corresponding 
to the three regions A, B and C as above, we have 

= ^/pa^ + y/qb^ + y/T^^J)^^c^ , where p, q and 1—p—q 
are the probabilities for an individual particle to be found 
in regions A, B and C respectively. We now trace the 
state 1^') = ^jy, |0) over the c-mode in order to find the 

reduced density matrix for modes d and b. Applying Eq. 
@) to the split AB — C it is clear that the reduced density 
matrix will correspond to the states 



(8) 



occuring with binomial probabilities (J^) [p 4- <?)"(! — p — 

First consider the limit where A and B are finite-size 
subsystems of an infinitely large BEC. The total parti- 
cle number N goes to infinity in such a way that the 
mean particle number in systems A and B together, 
A = N{p + q), remains finite. It is easy to show that in 
this limit the binomial distribution tends towards a Pois- 
son distribution. As a consequence the reduced state is 
exactly of the form (^) with the individual number states 
of the form (||). Thus in this limit there is no entangle- 
ment in the system consisting of A and B only. 

As long as C is of finite size, there is always entangle- 
ment between A and B. To show this, let us write the 
reduced state of modes d and b explicitly: 




T..i«)(*i=EE(:)a-!'-W(:) 

n=Ofc,i=0 ^ ^ V ^ ^ 

{V^)''+'{V^f''-''-%n- k){l,n - II (9) 
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where |fc, n — k) 



-|0). The inseparabiUty 



of (|^) can be shown using the Peres-Horodecki partial 
transposition criterion . 

This can be done by first projecting the state onto 
suitable two-dimensional subspaces in both modes. The 
simplest choice are the subspaces spanned by the states 
|0) and |1). The resulting state is 

(1 - p - <z)'^|00)(00| + 7V(1 - p - g)^-i(p|10)(10| 
+(Z|01)(01| + vw|io)(oi| + ^|01)(10|) 

+ ^^^^{l-p^qr-^2pq\ll){lll (10) 

whose partial transpose has a negative eigenvalue for all 
values of p and q, which is seen most easily by calculating 
its determinant. Thus any two regions within a finite-size 
BEC are entangled with each other. 

This is not the only possible choice of two-dimensional 
subspaces which exhibits entanglement. If one projects 
in both modes onto two arbitrary Fock states \n) and |m) 
(the same two numbers for modes d and 6), the resulting 
density matrix also has a non-positive partial transpose 
for all values of p and g, which can be shown in the same 
way as before. The amount of entanglement between the 
two regions A and B depends only on the total particle 
number N and on the average particle numbers in the two 
regions, determined by p and q. It does not depend on 
the distance between the two regions. This emphasizes 
the spatially coherent character of the BEC. 

We have noted above that the present type of entangle- 
ment cannot be detected by separate measurements on 
the modes d and h. However it can be detected by joint 
measurements of the two modes. A (conceptually) sim- 
ple possibility is to perform an interference experiment. 
For simplicity consider the state 



p = '^pMi){(l>i\® \xi){xi 
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[a) + b^) 



N 



-|o), 



(11) 



which has {ib\a) a 1^) = {ij}\b^b\i;) = N/2 and a fixed total 
particle number iV. Now consider detection in the new 
basis of modes given by a' = {a + b)/^/2, U = (d-6)/\/2. 
This corresponds to a measurement after superimposing 
modes d and 6 on a 50-50 beam splitter. The state ( |ll| ) 
satisfies {il)\a'^ a'\%j:) = N and {ipp^b']^^) =0. As a conse- 
quence of interference, behind the beam splitter all par- 
ticles are concentrated in one of the modes. 

There is no separable state of modes d and b which 
has the same properties, namely (i) & fixed total parti- 
cle number (different from zero), i.e. (AA^)^ = (iV^) — 
(7V)2 = 0, where N = a^a + b^b = d'^d' -I- and (li) 
complete destructive interference, i.e. {b'^b') — 0, at the 
same time. 

This can be seen in the following way. Consider a 
general separable state of modes d and b: 



where are states of mode d and \xi) are states of 
mode b. 

First assume that the particle number is fixed to be 
A^. This implies N\(j),)\xi) = N\(j)i)\xi) for all i. This is 
possible only if is an eigenstate of d^d and \xi) is an 
eigenstate of wb, i.e.|0i)|xi) = |ni)|A^ — m). So p must 
be a convex combination of products of Fock states. But 
Fock states have no definite phase, therefore destructive 
interference cannot occur. It is easy to see that in this 
case on has {b'^V) = TYpb'^V = N/2. The particles have 
a 50-50 distribution behind the beam splitter, in contrast 
to the result for (|ll|). 

On the other hand, one can also impose {b'%') = 0. 
This implies 

\Y.P, ((dtd), + (fotfe), - (at>,(6), - (d),(6t),) ^ 0, (13) 



where we have introduced the shorthand notation 
(atd)., =^ {(j)i\a^a\(j)i),{b^b)i = {xS'^b\Xi) ^ = 
((/)i|d|(/),), {b)i = {Xt\b\Xt) etc. 

Consider the expression in parentheses under the sum 
in (p^, dropping the index i for a moment. If we define 
(d) ^ ^e" and {b) = Be'l^ with A, B positive and a, [i 
real, this translates into 



A^ +Ca + B"^ + Cb- 2AB cos(a - (3). 



(14) 



Here we have introduced the correlation functions Ca = 
{a)d) — {a) {a)) and Cb = {b^b) — {b){¥), which are pos- 
itive definite, as can be seen from the Cauchy-Schwarz 
inequality apphed to vectors |0) and |?;) = d|0) |11|. 

Therefore (|lj) is never negative, and equal to zero only 
iiCa = Cb = 0,A = B and cos(a - /3) = 1. The Cauchy- 
Schwarz inequality is saturated only if the two vectors 
are collinear, therefore the first condition implies a\4>) = 
Ae"^\(p) and b\(f>) = Be''^\(t)). Returning to Eq. (|l|), this 
implies that all the states and \xi) have to be coher- 
ent states, with the additional constraints Ai = Bi and 
cos(ai — Pi) — 1 for all i. This implies that the parti- 
cle number cannot be fixed. For a mixture of products 
of coherent states of the form J2iPi\^i){^'i \ 
one has (TV) = (a^a + b^b) = 2j2,Pi\Ai\^ and (TV^) = 

(dtddtd+StS6tS+2dtd6t6) ^ 4.J2^pM^\^ + 2E^p^\M^ 

using dd^^ = d^d + 1. Since X^jPil^il* always greater 
than or equal to (X^jPil^iP)^' l^i^ implies that (AiV)^ = 
(A^) — (N)'^ > (TV). If one demands perfect destruc- 
tive interference, then (for a separable state) the particle 
number cannot be definite. 

Thus for separable states there is clearly a trade-off 
between the fiuctuations in particle number, expressed 
by (AA^)^, and the amount of interference, expressed by 
{b'^b'). If one of the two quantities is zero, the other one 
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is bounded away from zero. If these lower bounds are vi- 
olated, the state under consideration must be entangled. 
It may be possible to extend the above results, for exam- 
ple to derive a region in the ((AiV)^, {b'^b')) plane (for 
fixed mean particle number) where the states are defi- 
nitely entangled. This would be desirable because in a 
real experiment neither of the two quantities would be ex- 
actly equal to zero. Such a more general criterion might 
then allow a direct experimental proof of the presence of 
entanglement in Bose-Einstein condensates via interfer- 
ence. Note that interference between different distinct 
regions of a Bose-Einstein condensate was observed in 
|l2| , where the BEC was distributed over an array of 
micro-traps. 

The effect of non-zero temperature on the entangle- 
ment analyzed here is another topic for future research. 
More generally, one may wonder whether new insights 
in condensate phenomena such as superfluidity can be 
gained from the point of view of entanglement. Finally 
it should be interesting to compare the present results to 
the case of condensed fermion pairs. 

I would like to thank S. Bose and A. Costa for useful 
discussions. This work was supported by the QuComm 
project of the European Union (IST-1999-10033). 



ber superselection rule, which may be a fundamental law 
of nature or due to very fast environmentally induced 
decoherenc e. For a discussion of superselection rules cf. 
D. Giulini, 5uant-ph/0010090 Here these states are used 



only as a convenient mathematical tool for analyzing the 
properties of the mixed states (5) and (7). 
[9] Our present discussion is within the framework of non- 
relativistic quantum field theory. It is worth noting that 
in relativistic quantum field theory even the vacuum is 
an entangled state. See e.g. S.J. Summers and R. Werner, 
Phys. Lett. A 110, 257 (1985); L.J. Landau, Phys. Lett. 
A 123, 115 (1987); H. Halvorson and R. Clifton, J. Math 
Phys. 41, 1711 (2000); B. Reznik, |quant-ph/0008006 
(2000). 

[10] A. Peres, Phys. Rev. Lett. 77, 1413 (1996); M. Horodecki, 
P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1 
(1996). 

[11] {v\v} = {(p\a'' h\(t>) and {(t>\v) = {(p\a\(t>} fulfill the Cauchy- 
Schwarz inequality {v\v) {(j)\<l)} > {v\(j)) {(j)\v) . 

[12] B.P. Anderson and M. Kasevich, Science 282, 1686 
(1998). 



[1] J.S. Bell, Physics 1, 195 (1964); J.S. Bell, Speakable and 
Unspeakable tn Quantum Mechanics (Cambridge Univer- 
sity Press, Cambridge, 1987). 

[2] M.A. Nielsen and I.L. Chuang, Quantum computation 
and quantum information (Cambridge University Press, 
Cambridge, 2000). 

[3] D. Bouwmeester, A. Ekert, and A. Zeilinger, The Physics 
of Quantum Information (Springer, Berlin, 2000). 

[4] M.C. Arnesen, S. Bose, and V. Vedral, Phys. Rev. Lett. 
87, 017901 (2001); K.M. O'Connor and W.K. W ootters 
Phvs. Rev. A 63, 052302 (2001); W.K. Wootters, [quant 



ph/0001114| ; X. Wang, Phys. Rev. A 64, 012313 (2001); 
X. Wang, Phys. Lett. A 281, 101 (2001); D. Gunlycke, 
V.M. Kendon, V. Vedral, and S. Bose, Phys. Rev. A 64, 
042302 (2001). 
[5] C. Simon and J. Kempe, |quant-ph/0109102| (2001). 
[6] M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wie- 
man, and E.A. Cornell, Science 269, 198 (1995); K.B. 
Davis et al, Phys. Rev. Lett. 75, 3969 (1995); A.S. 
Parkins and D.F. Walls, Phys. Rep. 303, 1 (1998). 
[7] A. S0rensen, L.-M. Duan, J.I. Cirac, and P. ZoUer, Na- 
ture 409, 63 (2001) and K. Helmerson and L. You, Phys. 
Rev. Lett. 87, 170402 (2001) propose ways of creating 
entanglement between the internal degrees of freedom of 
particles in a BEC. 
[8] Such superposition states of different particle numbers 
are impossible to prepare in practice. If the particles are 
atoms or molecules, such states violate the baryon num- 



4 



